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Abstract. We calculate the heavy-to- light form factors in the relativistic quark 
model with heavy infrapropagators. Their g 2 -dependence in the physical region 
is defined by two parameters: the "infrared" parameter v characterizing the 
infrared behavior of the heavy quark and the mass difference of the heavy 
meson and heavy quark E = tuh — Mq. 

It is found that the values of the D — * K(K*) and D — ■> n(p) form factors 
at q 2 = are in excellent agreement with the available experimental data and 
other approaches whereas these values for B — * 7r(p) transitions are found to 
be larger than those of several other models. 

The obtained form factors are used to calculate the widths of the semilep- 
tonic decays of B and D mesons. The comparison of our results with the avail- 
able experimental data and other approaches is given. 



1 Introduction 



A theoretical study of the heavy-to-light transitions is more complicated than 
heavy-to-heavy ones because of the heavy quark symmetry cannot be applied in 
this case. The results of nonperturbative calculations made within lattice QCD 
[0-@: QCD sum rules and various quark models Q-|l3) are seemed to 

be rather different each other to give a reliable conclusion about the behavior 
of the form factors in the physical region. Nevertheless, such calculations may 
be useful for our understanding what is really important under description of 
the heavy-to- light transitions. 

One has to remark that the most of mentioned approaches allows one to 
calculate the values of the form factors at some fixed point of q 2 , the mo- 
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mentum transfer squared, and then to adopt some simple shape for getting q 2 
dependence in the full physical region. 

The lattice calculations give the form factors near q 2 = q 2 lax and then 
extrapolate them to q 2 = using polelike dependence. The QCD sum rules 
method allows one to give the reliable results at q 2 — and then also the 
pole dominance hypothesis is used for extrapolating them to physical region. 
The (^-dependence have been obtained by using the QCD sum rules method 
in the paper [ p7o[ . It was found that the pole dominance dependence is a good 
approximation for the vector form factors whereas this approximation is not 
applicable for the axial form factors. In the nonrelativistic constituent quark 
model fll| the form factors, computed with the overlap integral of the non- 
relativistic meson wave functions, have the exponential q 2 dependence. The 
quasipotential approach fl^| predicts similar behavior. The q 2 dependence of 
the relevant form factors was derived from the independent quark model fl3| 
with a scalar- vector-harmonic potential. Some nonrelativistic [jl4| and relativis- 
tic [t(| [TtJ quark models adopt polelike (monopole, dipole,...) ansatz for the 
q 2 dependence of the form factors. 

The dispersion relations incorporating with the Heavy Quark Effective The- 
ory (HQET) and Chiral Perturbation Theory (ChPT) ||, with the light-cone 
model [ [l9| , with the crossing symmetry and QCD perturbative calculation in 
unphysical kinematic region have been used for studying the form factors. 
In the paper pl| heavy-light form factors have been analyzed within approach 
incorporating HQET, ChPT and Vector Meson Dominance (VMD) model. 

Therefore it may be seen that a completely consistent calculation of the 
heavy-to-light form factors in the full physical region of the momentum transfer 
squared is a problem not only for the approaches based on the first principles 
of QCD but also for various quark models. The polelike q 2 dependence usually 
used for the form factors should be at least justified from the more deep physical 
representations. 

The present work is devoted to the calculation of the form factors describing 
the heavy-to-light transitions within a relativistic c 
light quarks and infrared heavy quark propagators 

Our previous approach |2^, Q was based upon the quark confinement model 
(QCM) which incorporates the confinement of light quarks by devising a quark 
propagator that has no singularities, thus forbidding the production of a free 
quark p^| . This model then allows one to perform covariant calculations of 
Feynman diagrams with dressed light quark propagators free of simple pole. It 
was shown that the QCM approach successful in many static and non-static 
properties of the light hadrons [Q . In the paper |g2| the infrared behavior of 
the heavy quark has been taken into account by modifying its conventional 
propagator in terms of a single parameter v. The weak decay constants and 
the Isgur-Wise function have been calculated. 

Here we slightly modify the inclusion of the infrared behavior of the heavy 
quark in a such way to provide the conservation of local properties of Feynman 
diagrams like gauge invariance and the Ward-Takachashi identity. We investi- 
gate the dependence of the heavy-to-light form factors on two parameters: the 



uark model with confined 
2l. 
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infrared parameter v and the mass difference of the heavy meson and heavy 
quark E = m# — Mq. 

We found that the values of the D — > K(K*) and D — > Tr(p) form factors 
at q 2 = are in excellent agreement with the available experimental data and 
other approaches whereas these values for B — > n(p) transitions are found to 
be larger those of several other models. 

The obtained form factors are used to calculate the widths of the semilep- 
tonic decays of B and D mesons. We compare our results with the available 
experimental data and other approaches. 



2 Model 



Our approach is based on the effective interaction Lagrangian which describes 
the transition of hadron into quarks. For example, the transition of heavy meson 
H into heavy Q and light q quarks may be described by 



£ in t{x)=g H H(x) dxx dx 2 6 



MQXi+m q x 2 



Mr 



f (xi-x 2 ) 2 Q(x 1 )r H \ H q(x 2 )+ii.c. 

(1) 

Here, Ah and /jj are the Gell-Mann and Dirac matrices, respectively, which 
provide the flavor and spin numbers of mesons H. The function f(x 2 ) charac- 
terizes the smearing of interaction between a hadron and quarks. 

The coupling constants gn defined by what is usually called the compos- 
iteness condition which means that the renormalization constant of the meson 
field is equal to zero: 



Z H = l-3g H /(2nffl' H (m 2 H ) 







(2) 



where n' H is the derivative of the meson mass operator. 

It is readily seen that in the heavy quark limit Mq ^> m q the Eq.(Q) becomes 

as 



Ant(a;) =>■ gHH{x)Q(x)r H X H dx 2 f 



(x - x 2 y 



q{x 2 ) + h.c. 



(3) 



which means the light degrees of freedom are factorized out from the heavy 
ones in according to the heavy quark symmetry. One has to remark that the 
function f(x 2 ) characterizing the long distance contributions is related to the 
light quark only in this limit. It leades to the modification of the light quark 
propagator in the momentum space: 



1 



m q - i> m q - 



(4) 



In the paper |2j| the monopole function f{p 2 ) = 1/(A 2 —p 2 ) has been used for 
calculations of physical values. But in this case the physical amplitudes had the 
threshold corresponding to quark production. In our approach pl[ an entire 
(nonpole) function has been used for a single light quark propagator to ensure 
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the quark confinement. 

1 f da-, 



m q — p j Az— p 1 

with the functions a and b being defined by 

2 f zda z 
J z l -p l 



a(- 



A 2 ' 



P u P ■ 



K-p 2 ) 



do. 



P 



(5) 



(6) 



To conserve the local properties of Feynman diagrams like the Ward identi- 
ties, we had used the following prescription for the modification of a line with 
n-light quarks within the Feyman diagram: 



n 



i 



m a - 



d<J z Yl 



1 



Az- 



-r, 



(7) 



In QCM we had used a simple choice of the confinement functions p4| 

a{u) = a,Q exp(— u 2 — aiu) b(u) = bo exp(— u 2 + b\u). (8) 

The parameters a.;, bi, and A have been determined from the best model de- 
scription of hadronic properties at low energies and the following values were 
found [Bl: 



ao = b = 2 



ai = 1 



h = 0.4 



and 



A = 460 MeV, 



which describe various physical observable quite well p4[ . 

Since a heavy quark (with mass Mq ) in a heavy meson is under the influence 
of soft gluons (which sets the scale Aqcd << Mq), it may be regarded as nearly 
on its mass-shell where the infrared regime should take place for its propagation. 
The infrared behavior for one-fermion Green's function (propagator) has been 
investigated in various papers (see, for instance, |2r| and the references therein). 
The result is well-known only for abelian gauge theories: 



S(p) ~ (to— py 



(9) 



where v ~ (a5/47r)(3 — A) with A being the gauge parameter. 

We have chosen [^2| the heavy quark propagator as in Eq.(||) to calculate 
the weak decay constants and the Isgur-Wise function. 

But as it is readily seen the using of such kind of propagators in the Feynman 
diagrams destroys some local relations like Ward identity which based on the 
obvious equality 
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1 



1 



1 



TO— ft- 



In this paper we suggest to take into account the infrared behavior of the 
heavy quarks by slightly different way which nevertheless allows one to conserve 



the local properties of the Feynman diagrams. We shall use the infrapropagator 
for a single heavy quark line in a diagram: 



{Mq- ft- f) 1+v ttv J u v du M Q +u- ft- i> 







drtsM+ it) (io) 

Here Mq is the constituent mass of a heavy quark and 

do* = -A»*?W*±± st(K +I f) = - l —— (11) 

ttv du uV ; MQ+u-ji-tf 

But we will use the prescription for a line with n heavy quarks analogous 
to those used in Eq.(Q) for averaging of light quarks: 

n Tjzmr^ =► k n u , r - ^ 



Let (Mq 1 , Mq^) be the structural integral corresponding to the 
Feynman diagram with n heavy quark local propagators. It is readily seen 



that the taking into account the infrared behavior according to Eqs.( [10| . [ll 12) 
gives 

K^(M Ql ,...,M Qn )= f d<j h u K^(M Ql +u,...,M Qn +u) 



sin(7r^) fdu d „tri\ nr „ , . sin(7rf) f r ^/ n \ .„ , „ . 

= — /— ^-K ( - HM Ql +u,...,M Qn +u)=^^\K ( -°\M Ql ,... > M Qn ) 

-kv ] U v OU TTV 





1 



-v 





■du K^{M Ql ,...,M Qn ) - K^(M Ql +u,-,M q „ + u) 



/^ (0 Ha%+«,-..,a^+")}-^^{^ (0) (Mq 1 ,...,m q j 

; du KW(M Ql , ...,M Qn ) - +«,..., M Q „ + u) 



i 

l 



(13) 



n 



This means that we can calculate the integral corresponding to Feynman di- 
agram with the local heavy quark propagators and then apply the averaging 
according to Eq.(O) to take into account the infrared behavior. 
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3 Heavy-to-light form factors 

First, we introduce the necessary notation. We will write P for light pseu- 
doscalar meson and V for light vector meson with the masses mp and my 
respectively. We will write H for heavy pseudoscalar meson with a mass mn ■ 

In what follows we assume that all masses and momenta in the structural 
integrals are given in units of A = 460 MeV. 

The invariant amplitudes for H — ► Iv and H — > P(V)lv transitions are 
defined as 

A(H - ev) = ^=V qq> (eO^)M£(p) (14) 
A(H -> Pev) = ^=V qq ,(eO^)M» HP {p,p') (15) 
A(H -> Ve^) = ^=V qq ,(eO^)e v M^ v {p,p'), (16) 

where 

M £W = ^9 H JdaZjda{J^M^ 1 5 S^+ p)0»S l z ()i)} = if H p» (17) 

= f + tfw+p'n + rtfw-p'n (is) 

M^(p,pO =^g H gvJdatJda l z J^tri { i 1 5 S^+ tf)0»S l z Ql+ p") 7 »S l z Q£)\ 

= -i(m H +m v )A 1 (q 2 )g^+i M{ - q ^ (jf+p'^+i ^ ) (p'-p'^p" 

mH+mv rriH+mv 

+ ^^ V ' (19) 
Here 

2 / l\2 2 2 /2 2 / nf ,\ 

q ={p-P) , P = ™> H , p =m P( y ) (20) 

The integrals like (0), (|l|), ([i||) are calculated by using the standard Feyman 
a - parametrization and the prescriptions (^) and ( p"3| ) to integrate over da l z 
and da^ , respectively. The explicit expressions for the form factors are given in 
Appendix. 

The hadron-quark coupling constants for light pseudoscalar and vector 
mesons and heavy pseudoscalar mesons determined from the compositeness 
condition (ffi) are written down 



2tt / 2 / 1 2tt / 1 

The functions Ri(x) are shown in Appendix. 



7 



4 Numerical results and discussions 

We have two adjustable parameters in our approach: the difference between 
the heavy meson and heavy quark masses E — — M c = raj — Mb and 
the parameter v characterizing the infrared behavior of heavy quark. We will 
adjust them by fitting the available experimental data for branching ratios of 
D-meson decays. The results of fit are given in Table 1. One has to remark 
that the allowed regions for the adjustable parameters are quite narrow: the 
parameter E varies from 0.315 up to 0.415 GeV whereas the infrared parameter 
v varies from 0.6 up 0.7. In what follows we will give our results for E=0.365 
GeV and i/=0.65. This set of the parameters seems to be best for describing 
the available experimental data. Our results for fo and fg in this case are 

f D = 193 MeV, f B = 136 MeV (22) 

The magnitudes of form factors of the D — ► K(K*)ev at q 2 — are given 
in Table 2. One can see that our results are in good agreement with the experi- 
mental data and other approaches. The dependence on the momentum transfer 
squared is shown in Fig.l and Fig. 2. 

The results for the form factors and decay widths of B-meson are given in 
Tables 3 and 4 by using the values of adjustable parameters obtained above. 
It is seen that the values of the form factors at q 2 = are substantially larger 
those obtained in QCD Sum Rules |fjjl- []lo| and quark models |^2[-|p!5| but close 
to the results of papers @, [§, §j|, |2J and ||| for /+(0). In the latter it was 
found that using the heavy quark symmetry and assuming pole dominance for 
the form factors, f+*(0) is estimated to be ss 0.39. If the requirement of heavy 
quark symmetry is weaken so that it applies only to soft pion emissions from 
the heavy meson, one finds /+"(()) ~ 0.53. 

Recently, the CLEO collaboration |29| claims the following results for 
semileptonic B decays branching ratios and \V u b\- 

Bv(B° -> tt-I+v) = (1.8 ± 0.4 ± 0.3 ± 0.2) x 10~ 4 

Br(S° -> p~l + v) = (2.5 ± 0.4t°;^ ± 0.5) x 10~ 4 (23) 

|K 6 | = (3.3 ±0.2^ ±0.7) x 10- 3 , 

One has to note that the branching fractions for B — > irev and B — > pev have 
been obtained by using a variety of models. What is given in Eq. (|2^) is just the 
average value of the branching fractions obtained for each model. The errors 
are statistical, systematic and the estimated model dependence. The values for 
|Vub| were extracted from the branching fractions using tb = 1-56 ± 0.05 ps. 
To obtain |Kib|aver, the tt and p modes were combined by fixing their ratio to 
the prediction for each model. 

This means that it is not easy to use these results for checking predictions 
of other models. Nevertheless, we extract the values r(B — > n(p)ev) /\V u b\ 2 x 
10 12 s _1 from Eq. ( p3| ) and put them without errors in Table 4. As one expects 
our predictions are almost twice larger those values obtained in a such manner. 



8 



At the same time the ratio ir/p is in a good agreement with the experimental 
result. 

The q 2 dependence of form factors is shown in Fig. 3. 

One can see from Fig. 1-3 that the behavior of the form factors f + (q 2 ) and 
V(q 2 ) in the physical regions of the momentum transfer squared {q 2 < 1.5 Gev 2 
for the D-meson decay and q 2 < 20.0 Gev 2 for the B-meson decay) agrees quite 
well with q 2 dependence of the monopole form factor. At the same time the 
behavior of the form factors Ai(q 2 ) and A2(q 2 ) is rather different from the 
monopole function and similar to the prediction of QCD sum rules Jic| . 

5 Conclusion 

We have slightly modified our previous work to take into consideration the 
infrared behavior of the heavy quark in consistent way to provide the conser- 
vation of local properties of Feynman diagrams like gauge invariance and the 
Ward-Takachashi identity. The infrared behavior of the heavy quark is mod- 
elled in terms of a single parameter v which modifies the simple free Feynman 
propagator. 

We have investigated the dependence of the heavy-to-light form factors 
calculated within this approach on two parameters: the infrared parameter v 
and the mass difference of the heavy meson and heavy quark E = ran — Mq. 

We found that the values of the D — ► K(K*) or D — > n(p) form factors 
at q 2 = are in excellent agreement with the available experimental data and 
other approaches whereas these values for B — > Tr(p) transitions are found to 
be larger than those of several other models. 

We have used the obtained form factors to calculate the widths of the 
semileptonic decays of B and D mesons. 
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Grant No. DE-FG05-84-ER40413, by the INTAS Grant 94-739 and by the 
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Two-point function: 



l 



fp 




o 



i 
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Rv{x) = B + -jdub(--^- > , 



Rh{x) 



do . J da 1 

(i 



Ua( ,)(o..r) ) + ( ^^L M 2 -±axjb(5(a,x) 



Heavy-to-light form factors: 



fn = 9H /r/fj;;/(/nn 



(l-«) 2 





+ ^(l-ai)(l-a2)m|f-M 2 + a2(l-ai)g 2 ^b^Zi(ai,a2)^ j 



Ma|Zi(ai,Q 2 ) 



/"(9 2 ) = -/ + (g 2 )+5«3p 



+ 



1 

Jdaa{l - 



(2^) 2 
M 2 



i i 

y dan |y y < 



cfa2 «2(1 — ai)m P b Zi(ai, 02) 



-4i(9 2 ) = 9agv 



(2n) 2 (m H +mv) 



da u < / daa 



M 



{I -af 



q b <5(a,5 ) 



+ /,/,,, /da 2 m2g+ ^~ g2 a(^(a 1; a 2 ) 




1 1 

+ Jdaijda2 





Ct\ 1— -Ql 



2 A(l-ai) 



m H — Qi 1 — 202+0102 



— a.2 I 1 — 011 +aia 2 — — a x a2 



1 ^(ai , a 2 )^ - Mb ^\(ai , a 2 )^ | 



M{q 2 ) = g H gv 



3 (m H + m v ) 
2(2tt) 2 



da u Jdaijda2 




1— 3Qi+2aia2+2af— 2af 02] a Z\ (01,02" 



+ Qi ^1 - 2ai - 2a 2 + 2aia 2 )^ Mb ^A(ai, a 2 )^j 
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1 1 

2ai(l - ai) a( A(on,oc2, 



A o(q 2 ) = -Ai (g 2 ) ~g H 9v ^-j0p^ / da^ / ' 1 ! / <>■ ■ :■. 



^(9 ) = gngv — 2 (2tt) 2 — / " J J 



l i 



ajMb Z\(ai,Q 2 ) + (l-a 1 )a[A(a 1 ,a 2 ) 



Here 

M — Mq + u 

oo oc 

A n = jdua(u)u n , B n — Jdub(u)u n , 8(a,x) = — ^M 1 — ax 

o o 

A(ai,a 2 ) = Ql — r-M 2 — ai(l — ai)m H — (1 — ai)(l — a2)a2mp (v) — QiQ2g 2 
(1 — Ql J 
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Table 1. Branching ratios for semileptonic and nonleptonic B and D meson decays 
for various values of the model parameters E and v. 



Decay 


Experimental 
average Pq| 


E = 0.415 
v = 0.60 


E = 0.365 
v = 0.65 


E = 0.315 
v = 0.70 


Br(D°^K-e+u e ) 


3.8 ±0.22 

fl oq+0.23 


4.26 
0.48 


3.77 
0.42 


3.32 
0.35 


Br{D°->K*-e+is e ) 


2.0 ±0.4 


1.39 


1.31 


1.23 


Br{D + ^K°e + v e ) 
Br{D+^iv l + vi) 


6.6 ±0.9 
0.57 ± 0.22 


11.0 
0.63 


9.82 
0.55 


8.57 
0.47 


Br{D + ->K*°e + v e ) 
Br{D + ^pe + v e ) 


4.8 ±0.5 
< 0.37 


3.54 
0.35 


3.34 
0.32 


3.14 
0.29 



Table 2. The predictions for the form factors at q 2 — in the decay D — > K{K* 
along with those of other approaches and the experiment. 



D^>K(K*)ev e 


/ + (0) 


Mo) 


.42(0) 


no) 


Expt. average 


0.75±0.04 


0.56±0.04 


0.4±0.08 


1.1±0.2 


Our results 


0.78 


0.52 


0.4 


1.07 



LMMS@ 

ELC|| 

APE§ 
UKQCD@ 
GSsf 
DP[| 




0-90^1 
0.63±0.0< l 

±15 
±7 



0.60 



0.78±0.08 

0.67 
0.71 
0.75±0.05 



+7 
8 

+12 +10 
12 -7 



0.83 ±2g 

0.53±0.03 
0.64±0.16 
0.67±0.11 
0.70j 



0.59^: 
0.19±0.21 



0.4 



±28 
= 4 



0.61 



!+7 

-10 
+6 +9 
-6 -7 



0.49±0.34 

0.66 
0.83 



+ 10 

15 

+20 +12 
20 -8 



0.86±0.1 
0.86±0.24 
1.08±0.22 
1.01±f 3 

1 -°^-24 -14 



C[2J| 



0.60±ig 


0.5±0.15 


0.6±0.15 


1.1±0.25 


0.76-0.82 


0.8 


0.8 


1.1 


0.73±0.07 


0.63±0.06 


0.43±0.04 


0.62±0.06 


0.8 


0.77 


1.48 


1.32 


0.7 


0.8 


0.6 


1.5 


0.76 


0.88 


1.15 


1.3 


0.78 


0.66 


0.43 


1.04 


0.67 


0.48 


0.27 


0.95 
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Table 3. The predictions for the form factors at q 2 — in the decay B — > ir(p)ev 
along with those of other approaches. 





/ + (0) 


MO) 


A 2 {0) 


V(0) 


Our results 


0.53 


0.50 


0.51 


0.70 
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Table 4. Semileptonic decay rates r(B — > n(p)ei/) in units 
comparison of our results with those of other approaches. 
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